On the ground states of the Bernasconi model 
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The ground states of the Bernasconi model are binary ±1 sequences of length N with low 
autocorrelations. We introduce the notion of perfect sequences, binary sequences with one- 
valued off-peak correlations of minimum amount. If they exist, they are ground states. Using 
results from the mathematical theory of cyclic difference sets, we specify all values of N 
for which perfect sequences do exist and how to construct them. For other values of N , 
we investigate almost perfect sequences, i.e. sequences with two-valued off-peak correlations 
of minimum amount. Numerical and analytical results support the conjecture that almost 
perfect sequences do exist for all values of N , but that they are not always ground states. 
We present a construction for low-energy configurations that works if N is the product of 
two odd primes. 



1 Introduction 

Binary sequences of +1 and —1 with low autocorrelations 
have many applications in communication engineering Q . 
Their construction has a long tradition , and has turned 
out to be a very hard mathematical problem. Bernasconi 
1^ introduced an Ising spin model that allows to formulate 
the construction problem in the framework of statistical 
mechanics. 

Consider a sequence of binary variables or Ising spins of 
length N, 



S = (so, si, . . . , sn~i) 
and their autocorrelations 

Af-l 



±1, 



(1) 



(2) 



where all indices are taken modulo N . Bernasconi defined 
a Hamiltonian H{S) by 

N-l 

H{S) ^ E^'(^) (3) 

9=1 

iV-1 Af-l 
t,j=0 g=l 

The ground states of this model with its long-range, four- 
spin interactions are the low autocorrelation binary se- 
quences we are looking for. 

The Bernasconi model is completely deterministic in the 
sense that there is no explicit or quenched disorder like in 



spin-glasses. Nevertheless the ground states are by defini- 
tion highly disordered. This self-induced disorder resem- 
bles the situation in real glasses. In fact, the Bernasconi 
model exhibits features of a glass transition like a jump in 
the specific heat ^ , and slow dynamics and aging Q . 

The replica scheme, an analytical method well estab- 
lished in spin-glass theory, is of little use for deterministic 
systems. People have approximated the Bernasconi model 
by a model with explicit, quenched disorder, which in 
turn can be analyzed analytically with the replica method 
1^, ^ . This approach is valid only in the high tempera- 
ture regime, however, and provides only estimates for the 
ground state energies. 

Bernasconi introduced his model in order to apply 
statistical mechanics, especially simulated annealing, to 
find low autocorrelation binary sequences. His approach 
turned out to be of little practical use, however: the en- 
ergy minima found with the annealing procedure differ by 
factor 2 from the conjectured ground state energy. This 
failure is due to the peculiar "golf course" property of the 
energy landscape which is dominated by a large number 
of local minima, while the global minima are deep and 
shallow holes, extremely isolated in configurations space. 

The original Bernasconi model is defined with aperiodic 
autocorrelations, i.e. with N ~ I ~ g as the upper sum- 
mation limit in equation (Q). For this model, exhaustive 
enumeration of all 2^ configurations seems to be the only 
means to get true ground states 

For the Bernasconi model with periodic autocorrela- 
tions, the situation is better. The tight connection be- 
tween the correlations of periodic binary sequences and 
mathematical objects called cyclic difference sets can be 
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exploited to get some exact results on the ground states. 
The theory of difference sets is well established in math- 
ematics but hardly known in statistical mechanics. It is 
one objective of this paper to fill this gap. It will turn out, 
that ground states of the Bernasconi model can be con- 
structed from cyclic difference sets for certain values of N. 
Other values of N require some generalizations. Guided 
by numerical results, we will discuss such generalizations. 

The paper is organized as follows: We start with a 
derivation of the equivalence of autocorrelations in binary 
sequences and cyclic difference sets in section ^ In section 
H, we introduce perfect sequences, which are ground states 
of the Bernasconi model - if they exist. The existence 
problem in turn can be answered using results from the 
theory of difference sets. Numerical and analytical inves- 
tigations of those cases, which are not covered by perfect 
sequences, can be found in section ^. Section ^ comprises 
our conclusions. 

2 Difference sets 

Let G = Zjv, the cyclic group of integers 0, . . . ,N— 1 
with addition modulo N, and let D C G be a /c-element 
subset of G. Each element g G G has a number of different 
representations g = di — 6,2 with di,d2 G D, denoted as 
replication number X{g): 

X{g) = \{{di,d2)\g^di-d2, di,d2£D}\. (4) 

The trivial case is A(0) — k. The values X{g) for g 7^ 
depend on D. The total number of differences di — d2 
equals fc^, leading to the constraint 



e = J2Xig). 

see 



(5) 



Now consider the periodic binary sequence (s^) associated 
with D: 



-1 i mod N e D 
-1 i mod N 1^ D 



(6) 



A straightforward calculation shows that the autocorrela- 
tions of this sequence are given by 



Cg^N-Aik-Xig)). 



(7) 



The converse is also true: Given a ±1 sequence of length 
A^ with periodic autocorrelations C{g), each g G G has 
X{g) = k + {C[g) — N)/A representations 5 = di — o?2 
within the fc-element set 



D:={i:Si = +l,i = Q,...,N-l}, 



(8) 



so we conclude that the periodic autocorrelations of a bi- 
nary sequence and the number of difference representa- 
tions of elements of a cyclic group G m D x D are equiv- 
alent. 



The definition of a replication number and its corre- 
spondence to the correlations of periodic binary sequences 
do not depend on the concrete realisation of our group G 
as G = Zat. Let G be a cyclic group of order A^, written 
in multiplicative notation. The replication number X{g) is 
defined as 



Mg) = |{('^i,c^2) I5 = di ■ d^^, di,d2 e D}\ 



(9) 



for a /c-element subset D C G. To construct a periodic 
binary sequence based on D, we choose a generator a of 
the cyclic group G and define the binary sequence to be 



+1 a' e D 
-1 a' ^ D ' 

Again the correlations of this sequence are given by eq. 
If the replication number is constant for g ^ 0, i.e. 



(10) 



Ha) 



k for 3 = 
A for 5 ^ 



(11) 



the set D is called an (A^, k, X) cyclic difference set, where 
the term cyclic refers to the fact that the underlying group 
G is cyclic - obvious variants are general or abelian differ- 
ence sets. D = {1, 2, 4} C Z7 is an example for a (7, 3, 1) 
cyclic difference set. 

Difference sets have been studied by mathematicians for 
more than half a century, see [p|-pT[ for a survey. They 
prove useful in various fields like finite geometry or design 
theory. According to eq. (|^) it is obvious that cyclic dif- 
ference sets and binary sequences with constant off-peak 
autocorrelations are essentially the same objects. 

The fundamental questions in the field of difference sets 
concern their existence and their construction. Albeit be- 
ing unsolved in general, the existence question can be an- 
swered for special sets of parameters {N,k,X). Most of 
the results are non-existence Theorems which impose nec- 
essary conditions on (A^, k, A) for a cyclic difference set 
to exist A simple non-existence Theorem is given by 
eq. (^), which for a (general) difference set reads 



fc(fc- 1) = X{N~1). 



(12) 



See |] for other, more sophisticated non-existence the- 
orems. All known existence proofs are constructive, i.e. 
they provide an explicit method to construct a cyclic dif- 
ference set. 

3 Perfect sequences 

We define a perfect sequence to be a binary sequence with 
one- valued off-peak autocorrelations of minimum amount, 
i.e. with 



Gg^Q - 



A^ = mod 4 

1 N = 1 mod 4 

2 A^ = 2 mod 4 
-1 N = 3 mod 4 



(13) 
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Eq. (^2|) excludes the existence of perfect sequences with 
all Ck>o = -2 for TV = 4t - 2 > 2. 

Obviously any perfect sequence is a ground state of the 
Bernasconi model. To find perfect sequences, we have to 
look for cyclic difference sets with parameters 

k~X^t = 4i,4t± l,4t- 2. (14) 

3.1 Hadamard difference sets 

Let us consider first the case = 3 mod 4. According 
to eq. (|l^), the cyclic difference set that corresponds to a 
perfect sequence must have parameters of the form 

(Af,fc,A) = (4i- l,2t- l,t- 1) (15) 

for some integer i > 0. Such difference sets are called 
Hadamard difference sets. All known Hadamard differ- 



N = 2^ — 1, j > 2 rn-sequence 

N = At + 3 prime Legendre sequence 

N = p(p + 2) p,p + 2 prime twin-prime sequence 

Table 1: Classes of A^ which allow the construction of a 
Hadamard difference set resp. a perfect binary sequence 
with all off-peak correlations Cg = —1. 



ence sets can be classified into three classes according to 
the value of A^ (table |^). In fact it has been shown that 
all Hadamard difference sets with A^ < 10000 belong to 
one of the above parameter classes, with seventeen possi- 
ble exceptions (all of them > 1000) jlj]. This result has 
been accomplished by the extensive use of non-existence 
Theorems. A physicist may conclude that there are no 
Hadamard difference sets for other parameters. The three 
classes cover 195 out of 250 values N ^ At + 3 < 1000. 
Table |^ lists a construction-rule for every parameter class. 
Beside these rules, there are additional construction meth- 
ods which lead to non-equivalent Hadamard difference sets 
for some parameter values like N ^ At + 3 ~ Ax^ + 27, A^ 
prime, for integer x 

We start with the description of m-sequences. Let p 
be a prime and F — GF{p) be the finite field of order p. 
A shift register sequence produced by a linear feedback 
shift register LFSR of order n over _F is a sequence a = 
(afc) of elements from F that satisfies the linear recurrence 
relation 

n 

Oi = (^jO-i-j i > n. (16) 

The sequence is uniquely determined by the initial 
conditions (ao,... ,a„_i) and the feedback coefficients 
(ci, . . . ,c„), Cj e F. The name "shift register sequence" 
stems from a hardware realization of eq. (p^). 



There are p" distinct n-tuples over F. The tupel 
(0, . . . ,0) reproduces itself under the linear recurrence re- 
lation. Therefore a LFSR of order n over GF{p) pro- 
duces a sequence which is ultimately periodic with period 
< — 1. A sequence with least period — 1 is called a 
maximum period sequence or, for short, an m-sequence. 

Note that in an m-sequence each n-tupel except occurs 
exactly once. Hence an m-sequence is independent of the 
initial conditions and solely determined by the recursion 
coefficients c. A well known Theorem ||l3|, ch.6] states 
that a LFSR with coefficients (ci, . . . , c„) produces an rn- 
sequence if and only if the associated feedback polynomial 

f{x) = 1 — CiX — C2X^ — ... — c„x" (17) 

is a primitive polynomial. Primitive polynomials of all 
degrees exist for finite fields of any order p, i.e. for every 
prime p there is an LFSR of length n that produces a 
sequence of period — 1, the maximum value. 

Supplied with an m-sequence over F — GF(2), we can 
define a binary sequence via 

S, = (-1)"' = (^_iyia.-l + -+c^a,-,^ _ (18) 

This sequence then satisfies the multiplicative recursion 

which produces an m-sequence, meaning that every n- 
tupel of ±l-values except the all-1 tupel (1, 1, . . . , 1) oc- 
curs exactly once as a subsequence (st, St+i, . . . , st+n-i) 
in (si). It is easy to see that the product sequence 
di — SiSi-^-g has the same recursion as the Si 

= dti ■ dT-2 rf?-n- (20) 

Since {si) is an m-sequence, {di) equals (s^) (for 5 7^ 0) 
except for a possible shift in the index, 

di = Si+s (21) 

for some s. Using this result, we find for the correlations 
of the sequence S = (sq, . . . , sn~i), A^ = 2" — 1, 

2" -2 2" -2 

since exactly 2"^^ of the Si equal —1 and the remaining 
2"-i - 1 equal -f 1. 

As we have seen, m-sequences over GF{2) with period 
2" — 1 can be used to construct perfect binary sequences 
with the same period. This construction requires a prim- 
itive polynomial of order n over GF{2). Tables of such 
polynomials can be found in the literature |p^-|l6| . 

Now we discuss the construction of Legendre sequences. 
Let p be an odd prime. Exactly half of the elements 
g e GF(p)* — GF{p) \ {0} are squares, i.e. they have a 
representation g = x'^ modp with x G GF(p)*. The other 
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half cannot be written as a square. A convenient notation 
for this property is the quadratic character 

r 5 = 

Xpig) = { 1 5 is a square in GF{p)* . 

[ —1 otherwise (23) 

The quadratic character can easily be calculated using Eu- 
ler's criterion chap. 15]: 



and it obeys 



Xp{g) =3^^^ modp, 



xeGF{p) 



(24) 



(25) 



for all 5 7^ 0. This orthogonality of Xp translates into good 
correlation properties of the so called Legendre sequence 
of length p 



1 i = mod p 
Xpi'i) otherwise 



(26) 



The periodic off-peak autocorrelations of the Legendre se- 
quence read 

p-i 

i=0 



Xp{g) + Xpi~g) + X! Xp{i)xpii + g) 
Xp{g) (i + (-i)(p-i)/2) - 1, 



„ , —1 p — i mod 4 
\ +1,-3 p=i mod 4 



(27) 



For p = 3 mod 4 the Legendre sequence corresponds to a 
cyclic Hadamard difference set. For p = 1 mod 4 exactly 
half of the autocorrelations take on the value +1, the other 
half -3. 

The discrete Fourier transform (DFT), 



N-l 



(28) 



3=0 



of a Legendre sequence is given by chap. 15] 



^ ( ^ - isk m^O mod p 
^ \ TO = mod p 



(29) 



for p = 3 mod 4. This peculiar similarity between the 
Legendre sequence and its DFT, 3mBk = — Sfe (fc ^ mod 



p), turns the Legendre sequences into ground states of 
another deterministic model with glassy properties, 



p-i 
fc=i 

p-i ^ p~i 

= ^ |sfc - — XI sm{27rjk/p)s 
k=i VP 



(30) 



2 



the so called sine-model 

The last entry in table |l| refers to the twin-prime se- 
quences. Let p, q be odd primes with p > q. Consider the 
binary sequence of length pq, 

+ 1 i = mod p 
— 1 i = mod q,i ^ mod p . 
Xp{i)Xq{i) otherwise (31) 

If p and q are twin primes [p — q -\- 2), all correlations of 
this sequence equal —1. A proof will be given in section ^, 
where we discuss the case of general odd prime numbers 
p and q. 

3.2 Menon difference sets 

Now we consider the case = mod 4. Eq. ( p^ ) enforces 
(TV, fc. A) = (4m^ 2ii2 - w, u2 - u) (32) 

for some integer w > for the cyclic difference set that 
corresponds to a perfect sequence. A cyclic difference set 
that satisfies eq. ( ^ ) is called Menon difference set. The 
corresponding binary sequence forms the row of a circulant 
Hadamard|^ matrix. Example for u ~ 1, written as binary 
sequence: 



(s) = (+1,+1,+1,-1). 



(33) 



Unfortunately, this is the only known example of a Menon 
difference set. In fact, it has been shown jlj] that no 
Menon difference set exists for 1 < u < 250 with the pos- 
sible exceptions u = 165 and u = 231. There is a tight 
connection between Menon difference sets and Barker se- 
quences, i.e. binary sequences with 



JV-l-g 

E 

i=0 



SiS 



i+g 



< 1 



(34) 



The existence of a Barker sequence of even length A'" 
implies the existence of a Menon difference set of order 
A^ |l^, It is not known whether the converse is also 
true. It has been shown, however, that there are no Barker 
sequences for 13 < A^ < 2.5 • 10^ with the possible excep- 
tion A^ = 1596961444 ||l8|. To summarize: It is a still 
unproven, but widely accepted conjecture, that no Menon 
difference set exists for A^ > 4. 



^This is the reason, why Menon difference sets are sometimes 
called Hadamard difference sets in the mathematical literature. Our 
ffadajuard difference sets are then called Paley-Hadamard difference 
sets p. 232]. 
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3.3 Other difference sets 

Another family of cyclic difference sets can be constructed 
for parameters 



(7V,fc,A) = 



q-l ' g - 1 ' q-l 



(35) 



where v < \ and g is a prime power. Such sets are called 
Singer difference sets [||. 

For q — 2we rediscover the parameters of the Hadamard 
difference sets based on m-sequences, q = 3,1^ = 1 corre- 
sponds to the lonely Menon difference set (^. The only 
other cases where Singer difference sets yield ground states 
of the Bernasconi model are 

q V {N,k,X) Cg^O 



4 1 (5,f,0) f 

5 1 (6,f,0) 2 
3 2 (13,4,f) 1 



Since the Singer difference sets yield ground states for only 
3 new values of N, we will not discuss their construction. 

According to eq. (p^, a difference set with N = I mod 4 
and Cg — 1 must have parameters of the form 

{N, k, A) = {2u{u + 1) + f , ^t^ - f )) (36) 



4.1 Almost perfect sequences 

An almost perfect sequence is a binary sequence with two- 
valued off-peak autocorrelations of minimum amount, i.e. 
with 



Cg-^o e 



{0,±4} iV = 0mod4 

{1,-3} iV = lmod4 

{-2,2} iV = 2mod4 

{-1,3} iV = 3mod4 



(38) 



Again we have an ambiguity: for = mod 4 all cor- 
relations are either G {0,-1-4} or e {0,-4}. A sequence 
with three-valued correlations Cg E {—4, 0, +4} does not 
match our definition. 

An almost perfect sequence corresponds to a set D with 
two-valued replication number where the two replication 
numbers differ by one: 



k .g = 
X{g) = A + 1 geU dG 
A giU 



(39) 



U can be any subset of the underlying cyclic group G. We 
call a set D with such replication number an almost cyclic 
difference set. 

Values of N that allow the construction of perfect se- 
quences are rare, as we have seen in the preceeding section. 
The question is whether the weaker constraint of almost 
perfectness can be fulfilled for more values of N . 



for some integer u > 0. It has been shown [gl| that no such 
cyclic difference set exists for 3 < u < 100, i.e. 13 < < 
20201, so again it is well-founded to believe that beyond 
TV = 13 there are no such sets at all. 

For N = At — 2, sequences with all Cfc>o = 2 require 



3t = -I- 2 



(37) 



with integer m, i.e. A^ — 2,6, 22, 26, 66, .... Except for 
N — 2 and A^ = 6, no such sequences are known. 

Other difference sets do exist, but either their param- 
eters coincide with one of the sets already discussed, or 
they lead to values for the off-peak autocorrelations that 
are far from a ground state. 



4 Beyond perfect sequences 

In the preceeding section we have seen that for many 
values of A^ perfect sequences do not exist. In this sec- 
tion we introduce a generalization called almost perfect 
sequences and discuss their appearance as groundstates of 
the Bernasconi model. In addition we present a construc- 
tion method for A^ = pq, p and q odd primes, that yields 
configurations with low energies, and we discuss general- 
ized Jacobi sequences. 



10° 



10' 



10 



10' 



10 



0.0 



□ □ 

A o 



A O 



'o o 



10.0 



O N=0 mod 4 
□ N=1 mod 4 
O N=2 mod 4 
A N=3 mod 4 



20.0 30.0 
N 



40.0 



Figure 1: Total number of almost perfect binary se- 
quences. 

To study this question we counted all almost perfect 
sequences of given length A^ by exhaustive enumeration 
(Figure |l|) Their number seems to increase exponentially 
with N but much slower than 2^. From figure |l| we may 
conclude that that there are almost perfect sequences for 
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all values of N but that their relative volume in configu- 
ration space decreases exponentially with N. 

For larger values of N, this conclusion is supported by 
a heuristic numerical search. We start with a random 
sequence and lower its energy by single spin-flips until we 
get a sequence whose energy can no longer be lowered by 
flipping a single spin. If this local minimum is an almost 
perfect sequence, the procedure stops. If not, it starts 
again with a new random configuration. 

With this algorithm we found almost perfect sequences 
for all N < 84. The CPU time increased strongly with N, 
so we stopped the search for larger values of N before an 
almost perfect sequence had been found. Random search 
is of course not well suited for configurations that are very 
rare. As an example, consider this: Nine workstations 
searched the configuration space for = 86 in parallel 
to find a sequence with Cg = ±2. After 20 days, two 
machines have found sequences with all Cg — ±2 for all 
g > except C43 — 6, while the other 7 machines still 
report sequences with at least two correlations Cg = 6 as 
their best results. The failure of random search for large 
values of N does not indicate the non-existence of almost 
perfect sequences. 

Compared to difference sets, much less is known about 
the construction of almost difference sets. We only know 
two construction methods for almost perfect sequences: 
Legendre sequences for N = At + 1 prime (see sec. ^) and 
a construction due to Lempel et al. |^ for A^ = — 1 
where p is an odd prime and to is a positive integer. 

The construction of Lempel et al. works like this: Con- 
sider the finite field F = GF{p"^). Let G denote the mul- 
tiplicative group of F and let a be any primitive element 
of F. The subset D of G defined by 



D = {a^ 



(40) 



where 



1) 



is an almost difference set (see [g2j for a proof) . The binary 
sequence that corresponds to D according to equation ( p^ 
is an almost perfect sequence with correlations 



Ca 



2 or —2 if k is odd 
or —4 if k is even 



(42) 



Example for p — 7 and to = 1: GF{7) is equivalent to 
Z7, the integers 0, . . . ,6 with addition and multiplication 
modulo 7. The primitive element a = 3 leads to 

^ = {2,5,4} (43) 

and the corresponding binary sequence reads 

5= (-1,-1, +1,-1, +!,+!) (44) 

with correlations Ci = C2 = C4 = C5 = — 2 and C3 = +2. 



For N = 2 mod 4, any almost perfect sequence is a 
ground state of the Bernasconi model, i.e. the construc- 
tion of Lempel et al. yields ground states for N = At + 2 = 
— 1, where p is an odd prime. Since our numerical 
results suggests that almost perfect sequences do exist 
for all values of A^, we arrive at the conjecture, that the 
ground states of the Bernasconi model for N ~ At + 2 are 
given by almost perfect sequences with ground state en- 
ergy E = 4(N — 1). This conjecture is in agreement with 
former numerical results j^]. 

For the three other residue classes, the question whether 
an almost perfect sequence is a ground state depends on 
\U\, i.e. on the number of correlations that deviate from 
the optimum value. 

Consider first the case A^ = 4i -I- 1 where i is a posi- 
tive integer. Let the value (— l)'(2i -I- 1) appear 2?!^ times 
among the correlations of a sequence, i — 0, 1, ... . The 
factor 2 refiects the fact, that Cg = Cn-q must hold. The 
possible excitations above the theoretical minimum energy 
E — At are given by 



AE 



1 1 °° 

iE-4t) = -J2^{^ + l) 



16 



— ni + 3^2 + 6n3 + 10^4 + 15?i5 



(45) 



Since X^i^o ^ must hold, the above sum is finite, of 
course. Equation (^ imposes the constraint 



C30 _^ 

^ j(7i2j - n2j-i) = -u{u + l) - t 



(46) 



on the rii, where u is a non- negative integer. An almost 
perfect sequence has Ui — for i > 1 and equation ( p6[ ) 
reduces to 



ni = t ~ 2^(^ + !)• 



(47) 



(41) This limits the possible values for ni to 



ni e {i,t-l,t-3, . . . ,i--Umax(Wmax + 1)} 



with 



c(t) = 



(48) 



(49) 



We ran a numerical search for low energy, almost perfect 
sequences up to size A^ = 81. Table |^ shows the result. For 
A^ < 41 we perform an exhaustive enumeration to check 
that the sequences we found are true ground states. Us- 
ing Eq. (^6|), we can assert that any hypothetical sequence 
with energy below the value listed in table || must be an 
almost perfect sequence for A^ = 45 and A^ = 49. This is 
not possible for larger values like A^ = 53, where a configu- 
ration with rii = and n2 = 1 has lower energy (Ai? = 3) 
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9 
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13 
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17 
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21 
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2 
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25 


3 5 6 


3 
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29 


14 6 7 


4 
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33 
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41 


4 7 9 10 


4 
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1 5 8 10 11 


(J 


Sb p 


49 


2 6 9 11 12 


6 


a p 


53 


3 7 10 12 13 


7 


7 


57 


4 8 11 13 14 


8 


7 


fil 


S 9 1 9 1 4 1 'i 


9 


7 




1 6 10 13 15 16 


1 


7 


69 


2,7,11,14,16,17 


11 


7 


73 


3,8,12,15,17,18 


12 


7 


77 


4,9,13,16,18,19 


18 


7 


81 


5,10,14,17,19,20 


19 


7 



quence must have 



Table 2: Number of correlations Cg = — 3 in low energy, 
almost perfect sequences for N = 4t+l found by numerical 
search. An entry a. p. in column "ground state" indicates 
that any hypothetical sequence with lower energy must be 
an almost perfect sequence. The energy of the sequences 
is given by E' = 16ni + iV — 1. 



and satisfies Eq. (|46|). For N — 5 and = 13, perfect se- 
quences based on Singer difference sets are ground states. 
Note that the Legendre sequences for = 1 mod 4 prime 
have ni = j{N — 1). This value is too large for a ground 
state, at least for those values of A^ considered in Table ^ 
Based on their numerical results for A^ < 50, Marinari 
et al. [D observed that 



AE = t-6 



(50) 



holds for i = 8, 9, . . . , 12. As can be seen from Table ||, 
this equation seems to hold at least up to i = 18 (A^ = 73). 

Consider now the case TV = 4t + 3. This time, 2ni 
denotes the number of correlations with value (— l)'"*"^ (2z+ 
1). The excitation energy reads 



1 1 

AE = _(£;-4i-2) = -^»(i + l)r 

i— 1 

= ?ii + 3^2 + 671,3 + 10n4 + 15n5 -f 
and Eq. (||) leads to the constraint 



(51) 



^ j(n2j - n2j_i) = -u(u + 1), 



(52) 



ni G {0, 1, 3, 6, . . . + 1)} (53) 



with 



(V9 + 16t- 1) 



(54) 



A 


ni: possible values 


found 


ground state 


27 


1,3,6,10 


3 


yes 


39 


1,3,6,10,15 


3 


yes 


51 


1,3,6,10,15,21 


6 


? 


55 


1,3,6,10,15,21 


10 


? 



Table 3: Number of correlations Cg = 3 in low energy, 
almost perfect sequences for A^ — At+3 found by numerical 
search. For all other values of A = 4t + 3 < 75, ground 
states are given by perfect sequences based on Hadamard 
difference sets. 

Most values A = 4t + 3 < 75 allow the construction of 
Hadamard cyclic difference sets, i.e. perfect sequences. Ta- 
ble ^ lists low energy sequences found by numerical search 
for the 4 remaining values N < 75. The numerical search 
for A^ = 75 yielded low energy configurations AE — 14 
with (711,712) — (8,2) and (ni,?i2) — (11,1). This does 
of course not exclude the possibility of an almost perfect 
sequence with lower energy. 

4.2 The case = mod 4 

So far all known ground states for A^ odd or A' = 2 mod 4 
are perfect or almost perfect sequences. The first excep- 
tions occur for A^ = At. Let rrii (rii) denote the number 
of correlations with value 4i (— 4i). The excitation above 
the theoretical minimum E — reads 



AE = ^E ^ V i^{m, + n,) 
16 

i—l 

and Eq. (||) leads to the constraint 

C30 

i{mi - n{) = u'^ -t, 



(55) 



(56) 



where u is a non-negative integer. An almost perfect se- 



where u is a non-negative integer. 

Table ^ shows low energy configurations in terms of rui 
and Ui, found either by exhaustive enumeration (A^ < 40) 
or numerical search. For A' = 28, 32, 36 and 40, the ground 
states are not given by almost perfect sequences, but by 
sequences which still obey \Cg\ < 4. For A^ — 48,56,64, 
even these bounds seem to be violated. 

The construction of Lempel et al. for A^ = p™ — 1, P 
odd prime, yields almost perfect sequences with ni = A/4. 
These sequences are not ground states for all values of A^ in 
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N 


A£' 


7712 


7Tll 


Til 


U9 
' 


ground state 


8 


1 








1 





VPS 


12 


1 





1 








ves 


16 


3 








3 





VPS 


20 


4 








4 





VPS 









4 








VPS 


24 


2 








2 





VPS 


28 


5 





1 


4 





VPS 








o 


2 


n 






4 




2 


2 


n 




40 


O 




2 


3 


n 




44 


Q 


n 


7 


2 


n 


? 


48 


6 


1 


2 








? 


52 


8 





2 


6 





? 




12 


I 


4 


4 


n 


? 


fiO 

\J\J 


7 


n 


4 


3 


n 


? 


64 


14 





6 


4 


1 


? 


68 


16 





12 


4 





? 


72 


13 





10 


3 





? 


76 


17 





7 


10 





? 



Table 4: Low energy configurations for N — At, found by 
exhaustive enumeration {N < 40) or numerical search. 



Table |. They do not even have the lowest energy among 
all almost perfect sequences. 

Wolfmann |2^ describes a construction for N = 
mod 4, that yields sequences with all Cg = 0, except 
Cn/2 = —{N — 4). The energy of such sequences is much 
higher than the ground state energy, at least for those N 
considered in Table ^ (except = 8). 

4.3 Two-prime sequences 

In this section, we describe the construction of low auto- 
correlation binary sequences for N — pq, p and q prime, 
which can be shown to be true ground states of the 
Bernasconi model if p and q are twin primes. For gen- 
eral primes p and q, they have at least a low energy, far 
below the energies that can be found by numerical search 
for large N. 

Let p and q be odd primes and p > q. Then G = 
Zq X Zp is a pq-element cyclic (additive) group. Consider 
the subset D C G given by 

D^ {ia,b)eG\a,b ^0, Xgja) ^ Xp(b)} 

' ' (57) 

U {(a,0) G G|a e ZJ 

D has 



k^-{p-\){q-\)-rq (58) 
elements and M is a multiplicative group. 



To calculate the correlations of the two-prime sequence 

+ 1 (i mod g, i mod p) ^ D , , 

— 1 otherwise ' 

we generalize the proof given in p5[ | to arbitrary primes 
p and q. Note that the sequences given by Eqs. (^) and 
(jsgj ) are the same. 

Let \{x, y) denote the number of different representa- 
tions 

(a;,?/) = (ai,5i) - (a2,62) {ai,bi) e D (60) 

of {x,y) e G — {0}. Since MD = D, differences {x,y) and 
(mix, 7712?/) with (ttii, 7712) G M will occur equally often in 
D. Hence there are constants Ai, A2, A3 and A4 such that 



A(x,y) 



Ai for x,y^O and Xq{x) = Xpiv) 

A2 ior x,y^O and Xgi^) Xpiv) 

A3 for a; = 
A4 



(61) 



for y = 

The corresponding autocorrelations read 

Cix,y)^4X{x,y)-pq + 2ip-q)-2. (62) 

To calculate Ai , . . . A4, we will look closely at the number 
of difference representations in Z^, the group of integers 
1, . . . , t— 1 with multiplication mod t, t prime. Let ei, £2 G 
{1,-1} and X e Zj and 

'^*i,e2(^) = \{{o-i^a2)\x = ai- a2. 



ai,a2 e Z*, Xi(ai) = £j|. (63) 



Then we have 



(64) 



We define A*^^^^ = ^\i,e2^^)- write -I-, — instead 

of 1,-1. Then clearly 



t- 1 



X\+ + X\_ + l = A*__ + A*_+ = 

From eq. (|2^) we deduce 

\\+ + \'__ - \\_ - A*_+ = -1. 
Furthermore, for e E {1, —1} we set 

et ^ / 1 if Xt(-l) = £ 
' ^ else 



(65) 



(66) 



(67) 



We decompose Zj in intervals of squares resp. nonsquares 
only, i.e. 



Zj = 1^ [rrii-i -I- 1, . . . , TTii], mo — 0, 771^ — t — I 



(68) 
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where [mi_i + 1, . . . , rui] consists of squares for odd i and and 
of nonsquares for even i. Then resp. A^j counts the 



difference representations 1 
i even, and thus 



(rrii + l) — mi for i odd resp. 



A*_+ - A^^_ = St 



(69) 



The above equations for the parameters A*^ provide an 
easily solvable system of linear equations with solution 



A* 



A* = 



1 

4 

t- 1 
4 

t~ 1 



SI 

Y 

s_*_ 
T 

SI 

T 



1 = 

t - 



t-3 



4 

t + 1 



(70) 



Finally, we define 



1 ifXp(-l) 
else 



1 if p — q mod 4 
else 



(71) 



For the value Ai we have: 



(72) 



Ai = A(l,l) 

= AV(A^_ + 1 + (5P) + A%(A% + J^) + 

Xl+iXl+ + 1 + ^^) + A1_(A^_ + ^^) + 1 + SP''^ 

where in the final term 1 + S^''^ the contribution 1 counts 
the difference representation (1,1) = (1,1) — (0,0) and 
SP-'J counts (1,1) = (0,0) - (-1,-1). Substituting the 
parameter values obtained above, we then obtain 

1 (5P^9 
X,^-{pq-2{p-q) + l) + —. (73) 

Let z e Z* with Xp{^) — ~1- We then obtain the value A2 
as: 

A2 = A(l,z) 

= X'L+{XP+_ + l + SP_) + Xl_{XP_++SP+)+ (74) 
A1_(A^+ + 1 + 6P^) + X%{XP__ +6P_) + 1- 6P'^ 

where the final term 1 — SP''' counts the difference rep- 
resentation (1,2:) = (0,0) — (— 1,— z). Substituting our 
parameter values gives 



A2 

Furthermore, 
A3 



1 SP''^ 
-{pq-2{p-q) + l)-—. 



(75) 



A(0,1) 

(A^++ + A^__) 



9-1 



9-1 



-(p+l)(g-l) 



(76) 



A4 = A(1,0) 
V9 



(AU + A1_)- 



2 +X\_+X% + 2 



i(g-3)(p-l) + g. 



(77) 



The correlations for 5 7^ finally read 

1 + 2(5P'« for Xq{g) = Xpig), .9^0 mod p, q 
C={ 1-2(5P'« for Xq(.9) 7^ Xp(5) 5 ^ modp,(7 

(78) 



^ ^ p — q — 3 for 5 = mod g 
q — p + 1 for (? = mod p 



with 



£;^(l + 45P'«)(p- 1)((7- 1) + 

(p-g-3)2(p-l) + (p-9-l)2(g-l). 



(79) 



For p — q = 2, all correlations are —1, i.e. we proved that 
the twin-prime sequence Eq. ( ^ ) corresponds to a cyclic 
Hadamard difference set. For p — g > 2, the construction 
does not necessarily lead to a ground state, as can be 
seen for N = 21, 33, 57, 65, 69 (table |) and N = 39, 51, 55 
(table ^) . Nevertheless two-prime sequences have energies 
well below those found by extensive numerical search for 
large N = pq and small difference p — q. 

4.4 Jacobi sequences 

The Jacobi symbol iPnU) is a generalization of the 
quadratic character XpU) (Eq- ^) to the case that N ~ 
ViV'i ■ ■ 'Ps" is the product of odd primes: 



V'jv(j)=n4".w- 



(80) 



Note that i/'w(i) = if (j,_Pi) = 0, i.e. if j is divided by 
Pi . The number of such zeroes in the range ... A'^ — 1 is 
given by — (piN), where Euler's totient function (p{N) 
is defined as the number of positive integers smaller than 
N that are coprime to N. In our case we have 



ct>{N)='[[pr\p^-l) 



(81) 



i=l 



A Jacobi sequence is a ±1 sequence of length N, given by 



iPnU) 
±1 



if VAr(j) y^O 
if V'Ar(j) =0 



(82) 



Setting the zeroes of the Jacobi symbol to ±1 leads to 
2N-<t>{N) (;jistinct sequences. 

Jacobi sequences are potential candidates for low en- 
ergy configurations, for three reasons: First, the Jacobi 
sequences reduce to the Legendre sequences for N prime. 
Second, for N = piP2 and a special choice for the zeroes of 
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iN 


P 


q 


two-prime 


Energies 
Jacobi 


ground state 


21 


7 


3 


84 


84 


42 


33 


11 


3 


228 


208 


64 


39 


13 


3 


774 


342 


86 


51 


17 


3 


2306 


610 


< 146 


55 


11 


5 


230 


230 


< 214 


57 


19 


3 


3672 


648 


< 184 


65 


13 


5 


736 


592 


< 224 


69 


23 


3 


7300 


964 


< 244 


77 


11 


7 


364 


364 


< 364 


85 


17 


5 


2100 


1076 


< 324 


87 






1 61 1 8 


1 71 8 


< 438 


91 


13 


7 


330 


330 


< 330 


93 


31 


3 


20508 


1788 


< 764 


95 


19 


5 


2926 


1102 


< 894 



Table 5: Energies of configurations for N ~ p*q, p,q odd 
primes and p > q + 2. For p — q-\-2, twin prime sequences 
are ground states. The column "Jacobi" contains the min- 
imum enery of all Jacobi sequences for given iV, obtained 
by exhaustive enumeration of all 2^~'^(^) = 2^'+''"^ such 
sequences. 



the Jacobi symbol, we rediscover the two-prime sequences, 
eq. (^). Third, Borsari et al j2^] have shown, that the 
Jacobi sequences are ground states of the sine model, eq. 
(|30|), for = 4t -I- 3 and not divisible by a square. 

To check the suitability of the Jacobi sequences as low 
energy configurations for the Bernasconi model, we did an 
exhaustive enumeration of all 2^~'^(^) Jacobi sequences 
of fixed length N = pq to find the sequence with mini- 
mum energy. Table |^ displays the result. The two-prime 
sequences are special Jacobi sequences, so their energy can 
not be lower than the Jacobi energy. The data from table 
1^ indicates, that the two-prime sequences have minimum 
energy among all Jacobi sequences for p—q < 6. For larger 
values of p ~ q, the energy of the two-prime sequences is 
larger. On the other hand, the discrepancy between the 
Jacobi energy and the ground state energy is clearly visi- 
ble. Except for the twin-prime case, Jacobi sequences are 
not ground states of the Bernasconi model. This discrep- 
ancy is even larger for non-squarefree values of N (Table 

I). 

This is no surprise, since prime powers dividing TV lead 
to regularities in the Jacobi symbol which keep the corre- 
lation values high. As an example consider N — p^: In 
this case is periodic, ipNij + p) = i/^nU), and we get 

Ck > p'''\p - 2) if A: = mod p (83) 

for all Jacobi sequences. 



AT 

IN 


T 1 • 

Jacobi 


ground state 


9 


56 


24 


25 


3000 


72 


27 


3802 


74 






< 1 94 


49 


2055 


< 144 


63 


2622 


62 


75 


17354 


< 298 



Table 6: Minimum energies of Jacobi sequences for odd, 
non-squarefree values of iV < 100, compared to the ground 
state energies of the Bernasconi model. 



5 Summary and open problems 

We have demonstrated, that the ground states of the 
Bernasconi model with periodic boundary conditions are 
closely related to cyclic or almost cyclic difference sets. Us- 
ing theorems and well founded conjectures from the theory 
of cyclic difference sets together with exact enumerations 
and numerical searches, we get some results on the ground 
states of the Bernasconi model, which are conveniently 
summarized by the following facts and conjectures. 

Fact 1 Perfect binary sequences, Eq. ( f?^ do exist and 
form ground states of the Bernasconi model for N = 
4,5,6,13, iV = 2-'' - 1 with j > I, N ^ At + S prime, 
N — p{p + 2) with both p and p + 2 prime. 

This fact can be proven using the known constructions for 
cyclic difference sets, as has been shown in section |^. 

Conjecture 1 For all values of N that do not match any 
of the conditions in Fact there is no perfect sequence, 
i.e. ground states are formed by non-perfect sequences. 

This conjecture has a very solid base in the various non- 
existence theorems for cyclic difference sets as cited in sec- 
tion |. 

Conjecture 2 Almost perfect sequences exist for all val- 
ues of N . 

Beside those values of N that match a condition in Fact 
|ll this can be proven for N ^ At + 1 prime (Legendre 
sequence) and N ^ q — 1, where q is an odd prime power 
(construction of Lempel et al.). For other values of N, we 
checked this conjecture numerically up to iV < 84. 

Conjecture 3 For N = 2 mod 4, the ground states of the 
Bernasconi model are given by almost perfect sequences. 
The ground state energy is E = 4(iV— 1). For other values 
of N , the ground state energy obeys 

( 16(iV-l) /or iV = mod 4 
-\9(iV-l) for N = 3 mod 4 ^ ' 
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This follows immediately from conjecture ^. 

Conjecture 4 For N = 1 mod 4, almost perfect se- 
quences are ground states of the Bernasconi model. The 
ground state energy is E — 5{N — 1) — 8u(u + l), where u 
is a non-negative integer < ^\/2N — 1. 

This conjecture is based on the numerical results from 
Table For the thermodynamic Hmit, our results lead us 
to 

Conjecture 5 The ground state energy of the Bernasconi 
model is 0{N). The limit limjv->oo "^/^ "is however not 
well defined, since the precise value of E{N) depends on 
number theoretic properties of N . 

This result is very different from the situation for the 
Bernasconi model with aperiodic boundary conditions, 
where E{N) = 0{N^), and limAr^oo ^pJP- seems to be 
well defined 0. 

Marinari et al. Q noticed some patterns in their numer- 
ical results for the ground state energies for < 50. Their 
equation (12) is a special case of conjecture |[ Eq. (13) of 
ref. ^ is equivalent to conjecture ^ with fixed value u = 3 
for N > 33. Our numerical results (Table |^) support this 
special value of u up to iV = 73. For iV = 77 and iV = 88 
we only found u < 1- 

We did not find sequences with lower energies than those 
reported in ref. This is due to the fact that numeri- 
cal methods like the one applied by Marinari et al. work 
quite well if N is not too large {N < 50). Our experience 
with numerical search suggests, that beyond N > 100 true 
ground-states can only be found by mathematical insight 
rather than computer power. 

One step towards more mathematical insight is to prove 
or disprove the above conjectures. Especially the conjec- 
ture on the existence of almost perfect sequences for all 
values of N deserves some attention. Another open prob- 
lem is the generalization of successful construction meth- 
ods for special N to more general values. Here, the method 
of Lempel et al. ||2^ is the most interesting candidate. 

Acknowledgement: The authors appreciate fruitful 
discussions with Alexander Pott. 
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